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ly-^ ■ The Bardeen-Cooper-Schrieffer (BCS) mean-field theory of the pairing interaction breaks 
down for nuclei and ultra-small metallic grains (nanoparticles). Finite-temperature pair- 
ing correlations in such finite-size systems can be calculated beyond the BCS theory in 
an auxiliary-field Monte Carlo approach. We identify thermal signatures of pairing cor- 
1^ relations in both nuclei and nanoparticles that depend on the particle-number parity. 
O ' 

g ; 1. INTRODUCTION 

Q . The pairing interaction in bulk metals leads to a phase transition from a normal to a 
superconducting metal below a certain critical temperature. This phase transition was 
explained by the Bardeen-Cooper-Schrieffer (BCS) mean-field theory yy. BCS is valid 
(J • when the pairing gap A is much larger than the single-particle mean-level spacing S, a 
3 ■ condition that is satisfied in bulk metals. 

The pairing interaction also plays an important role in nuclei. Pairing effects in nuclei 
at zero temperature are well documented but much less is known about thermal signa- 
^ I tures of pairing correlations. In nuclei A/ 5 is of order 1, and a typical nucleus is in the 
' crossover regime between the bulk BCS limit and the fluctuation-dominated limit. In this 
crossover regime, fluctuations tend to wash out signatures of the pairing transition and an 
interesting question is whether pairing correlations survive despite the large fluctuations. 

Experimental and theoretical studies of ultra-small metallic particles (nanoparticles) 
have shed light on pairing correlations in flnite-size systems These ultra-small metallic 
grains can be connected to leads and their transport properties measured. The number 
of electrons on the grain is controlled by changing a gate voltage. Nanoparticles whose 
linear size is below ~ 5 nm are also close to the fluctuations-dominated regime. 

We have studied pairing correlations in nuclei beyond the BCS theory using auxiliary- 
fleld Monte Carlo (AFMC) methods in the framework of the nuclear shell model. These 
techniques are also known as the shell model Monte Carlo (SMMC) methods [3.4j. With 
AFMC we can carry out fully correlated and realistic calculations in much larger conflgu- 
ration spaces than those that can be treated by conventional methods. Recently, we have 
extended the AFMC methods to study pairing correlations in nanoparticles beyond the 
BCS approximation [5]. 

The AFMC methods are briefly discussed in Sec. |2l Signatures of pairing correlations 
in the heat capacity piT] and moment of inertia of nuclei [8], calculated in the AFMC 
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approach, are presented in Sec. [31 Analogous signatures of pairing correlations in the heat 
capacity and spin susceptibility of nanoparticles [5] are presented in Sec. ID 

2. AUXILIARY-FIELD MONTE CARLO (AFMC) METHODS 

Correlations beyond the mean-field approximation can be calculated by taking into 
account fluctuations of the mean field. This can be formally expressed by the Hubbard- 
Stratonovich (HS) transformation [9]. The Gibbs ensemble e~^^ at inverse temperature 
(5 = 1/T (which is also the many-body propagator in imaginary time (3) can be written 
as a functional integral 

e-P^ = f V[a]GMa (1) 



over propagators describing non-interacting particles in time-dependent fields o"(r). 
The quantity in ([1]) is a Gaussian weight. 

Finite-size effects are important in both nuclei and nanoparticles. For such small sys- 
tems, it is necessary to calculate thermal expectation values in the canonical ensemble in 
which the number of particles A is fixed. For an observable O 

where Tr^ denotes a trace at a fixed number of particles A. = G^\TtaU(j\ is a 

positive-definite weight function, = TrA?7cr/|Tryi?7o-| is the Monte carlo "sign" and 
{0)a = TTA[OUo-]/TTAUa- Both Ti aU^ and {0)a can be evaluated using matrix algebra 
in the single-particle space. For example, the grand canonical trace of the one-body 
propagator U„ is given by 

TtU^ = det(l + U,) , (3) 

where Uo- is the A^^ x Ng matrix representing 11^ in the single-particle space containing 
Ng single-particle orbitals. 

The canonical trace can be evaluated using particle-number projection with 0^ = 
2Tmi/Ns as quadrature points. For example, the canonical trace of Uc, is given by 

Tr^f/. = E e"''^'"^ det (l + e'<^-U.) . (4) 

The multi-dimensional integral over the auxiliary fields in ([1]) is evaluated by Monte 
Carlo methods. The auxiliary fields are sampled according to the distribution W^- For 
samples {cXi} the expectation value in ([2]) is estimated from 

(O) ^ ^^1^^ . (5) 

Such auxiliary field Monte Carlo (AFMC) methods have been developed in the frame- 
work of the interacting nuclear shell model [lOj. We have recently extended AFMC to 
nanoparticles [5]. 
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For a general interaction, the sign can fluctuate from sample to sample. When 
the statistical error of the sign is larger than its average value, the method breaks down. 
This is known as the Monte carlo sign problem. In the nuclear case, the dominating 
collective components of the interaction have a good sign. Interactions with small bad- 
sign components can be treated by the method of Ref . [1] . 

3. NUCLEI 

We have used AFMC to calculate statistical properties of nuclei in the iron region 
within the complete fpgQ/2 shell [TU]. The single-particle Hamiltonian corresponds to 
a Woods-Saxon potential V plus a spin-orbit interaction [11]. The two-body interac- 
tion [T2IT3] includes a monopole pairing interaction whose coupling strength is determined 
from the experimental odd-even mass differences. Also included are multipole-multipole 
interactions that are obtained by expanding the surface-peaked interaction v{r, r') = 
—x{dV/dr){dV/dr')6{r — r') into quadrupole, octupole and hexadecupole components. 
The coupling constant x is determined self-consistently [12] = /o°° dr (dV/dr) (dp/dr) 
(p is the nuclear density). The quadrupole, octupole and hexadecupole interactions are 
then renormalized by 2, 1.5 and 1, respectively. All interaction components are attractive 
and lead to a good-sign Hamiltonian. 

In the following subsections we discuss signatures of pairing correlations in the heat 
capacity and moment of inertia of nuclei. 

3.1. Heat capacity 

The heat capacity C = dE/dT is calculated by a numerical derivative of the thermal 
energy. We used the method of Ref. p|, in which the thermal energy at f3 ± 6f3 are 
calculated with the same Monte Carlo walk and correlated errors are taken into account 
to estimate the statistical error of the numerical derivative. 

The Monte Carlo approach in a single major shell is valid for temperatures up to 
T ~ 1.5—2 MeV. At higher temperatures the heat capacity saturates because of truncation 
effects and it is necessary to include higher shells. Monte carlo calculations in larger spaces 
are possible but may be time consuming. Instead, we have combined the fully correlated 
calculations in the truncated space with independent-particle model calculations in the 
full space (including all bound states and the continuum) [7]. 

The right panels of Fig. [1] show the heat capacity versus temperature in the even- 
even nucleus ^°Fe and even-odd nucleus ^^Fe. The solid squares describe the extended 
heat capacity, calculated by the method discussed above. These heat capacities have the 
correct approximate linear behavior at higher temperatures. The dashed lines in Fig. [T] 
are the results of the independent-particle model. In the even-even nucleus we observe 
an enhancement in the heat capacity (an 5* shape curve), while in the even-odd nucleus 
the heat capacity remains close to the heat capacity of the independent-particle model. 
Similar even-odd effects were measured in the heat capacity of rare-earth nuclei [1A\ as is 
demonstrated in the left panels of Fig. [TJ 

3.2. Moment of inertia 

The moment of inertia / at finite temperature determines the spin distribution of energy 
levels in the framework of the spin cut-off model. It describes the response of the nucleus 
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Figure 1. Right: Heat capacities of ^°Fe (top) and ^^Fe (bottom) versus temperature T [7]. 
The sohd squares describe the extended heat capacities while the dashed hues describe 
the heat capacities in the independent-particle model. Left: measured heat capacities 
in ^^^Yb and ^^^Yb (solid lines) are compared with the heat capacities in the Fermi gas 
model (dashed lines) [14j . 



to rotations. For a rotationally-invariant Hamiltonian, / is given by / = ), where Jz 
is the 2;-projection of total angular momentum. 

Figure [2] shows the SMMC moment of inertia (symbols) as a function of f3 for a series 
of iron isotopes. In general, we observe a suppression of the moment of inertia at low 
temperatures because of pairing correlations. In the even-odd isotopes, this suppression 
is weaker and we even observe an enhancement in the limit T — > 0. 

This behavior of the moment of inertia can be explained in a simple model of a de- 
formed Woods-Saxon potential plus a monopole pairing interaction [S]. The single-particle 
states are organized in time-reversed pairs k and k (corresponding to angular momentum 
projection m and — m along the symmetry axis). The simple Hamiltonian is then given 

by 

-f^def = (^k{alak + a|afc) - gP^P , (6) 

k>0 

where = J2k>o '^I'^s is the pair creation operator and g is the pairing strength. 
The major odd-even effects are captured by a number-parity projection 



P, = 1(1 + Ve'""^) , 



(7) 
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Figure 2. The moment of inertia versus T for a family of iron isotopes. The circles are the 
SMMC results, while the lines are the results of a simple model (see text). The solid lines 
are the even number-parity projection while the dashed lines describe the odd number- 
parity projection (for neutrons). The dotted-dashed lines are the rigid-body moment of 
inertia. The deformation parameter is (32 = 0.14 and the pairing strengths are Qp = 0.42 
and Qn = 0.36 for protons and neutrons, respectively. From Ref. [8]. 



where A is the particle-number operator. projects on even (odd) number of particles 
for = 1 (r/ = — 1). 

For a deformed nucleus, rotational symmetry is broken and the moment of inertia is a 
tensor lij. The number-parity projected intrinsic moment of inertia is given by 

^ ^ lo^dT (J,(r)J,(0)) + r/nfc>otanh^f"/o^rfr ( J.(r) J,(0)). 

l + r^n.>otanh^^ ' 

where Jiir) describes the i-th component of the angular momentum in the intrinsic frame 
at imaginary time r. The expression for Jq dr {Ji{T)Jj{0))T^ is obtained from the expression 
for Jq dr {Ji{T)Jj{0)) (which depends on the quasi-particle occupations fk = jiji^b^) by 
the substitution 

fk^fk= i_^pE, ■ (9) 

Fluctuations in the pairing order parameter (i.e., the gap) are included in the static 
path approximation (SPA). Finally, rotational symmetry is restored by integrating over 
all orientations of the intrinsic frame. The results of the simple model are shown by the 
lines in Fig. [2l 
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4. NANOPARTICLES 

The spectra of ultra-small metallic grains (nanoparticles) of size ~ 2 — 10 nm were 
determined in Refs. |15J by connecting them to leads and measuring their non-linear 
conductance. The number of electrons on the grain is controlled by varying a gate voltage. 

The Hamiltonian of the metallic grain is given by 

-f^nanoparticle = H ( " eaVg + -Qb^-bBg ) fl^'^Aa + ^TT^ 4+'^A-'^At-'^M+ • (^O) 

A<7 ^ / ZL. 

The first term on the r.h.s. is a one-body Hamiltonian describing electrons in spin- 
degenerate orbital levels A (a = ± corresponds to spin up/down electrons). Also included 
is a Zeeman term describing the coupling to an external magnetic field B (/xb is the Bohr 
magneton and (^s is the (yf-factor). e^iV^/2C is a charging energy term that describes the 
Coulomb energy of the grain with capacitance C. The last term on the r.h.s. is a reduced 
BCS pairing interaction describing the scattering spin up/down electron pairs between 
different orbital levels. 

The Hamiltonian flTU]) can be solved in the BCS approximation. Below a certain critical 
temperature T^, a superconducting solution exists. BCS theory is a mean-field approx- 
imation and is valid in the limit when the pairing gap A is large compared with the 
single-particle mean-level spacing 5, .i.e., A/(5 ^ 1. This holds for the larger grains and 
a pairing gap was observed in the excitation spectrum of a grain with an even number of 
electrons P^ . 

However, for the smaller grains the ratio A/ 5 is comparable or smaller than 1. This is 
the crossover to the fluctuation-dominated regime where BCS theory breaks down. Static 
fluctuations in the gap order parameter can be taken onto account in the SPA and small 
amplitude quantal fluctuations are included in the SPA plus RPA approach |T6]. However, 
at lower temperatures it is necessary to account for additional quantal fluctuations. For 
small number of single-particle levels and electrons, diagonalization methods were used 
to find the eigenvalues [17]. For the pairing Hamiltonian, it is possible to reduce the 
eigenvalue problem to a set of non-linear equations [18], and this method was used in 
Ref. [19]. However, at higher temperatures the number of relevant many-body levels 
increases rapidly and the method becomes impractical. 

We have recently extended the AFMC methods to nanoparticles [S]. The attractive 
pairing interaction is a good-sign interaction in a density decomposition and thus accurate 
Monte Carlo calculations are possible. The AFMC method can be used for a large number 
of single-particle levels and are valid at both low and high temperatures. A different Monte 
Carlo method that is suitable for pairing-type interactions was also applied recently to 
nanoparticles [20] . 

4.1. AFMC for nanoparticles 

In the following we consider a grain with a fixed number of electrons in the absence 
of magnetic field. The charging energy becomes a constant and can be ignored. The 
Hamiltonian ffTOj) can be written in a density decomposition 



H = Y, ^xfix - - ^[(pam + Pam)^ - (Pam - Pa^.)^] , 

A A/i 



(11) 
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where n\ = a\^ax^ + a\_ax- is the number operator of level A, pxfj, = a\_^_a^+ are spin-up 
density operators and px^ = at^_a^_ is the time-reversed operator of pxfj.- The interaction 
in ffTTj) has a quadratic form and is thus suitable for an HS decomposition ([1]). 

In a metallic grain, the band width is determined by the Debye frequency ujo- In 
an Al grain, the measured Debye frequency is ujd ~ 34 meV. For the gap, we use the 
experimental value of thin films A ^ 0.38 meV. The coupling constant is determined from 
A ^ 2ujDe-^'^ to he 0.193. In a grain with A/(5 = 1, the number of single-particle 

levels is then k, 180. 

In practical AFMC calculations, we truncate the model space to a smaller band width. 
This can be done by renormalizing the coupling strength g. A simple way of determining 
the new coupling constant in the truncated space is by keeping the BCS gap (for a 
discrete single-particle spectrum) fixed. For a picketfence (i.e., equally spaced) spectrum, 
the renormalized constant Qr in a truncated space with 2Nr + 1 single-particle levels and 
half filling is given approximately by 



The choice of the truncated model space depends on the temperature range of interest 
and higher temperatures require larger model spaces. A model space of = 25 is good 
for temperatures up to T ~ 3S. The AFMC calculations scales as ~ and at lower 
temperatures we can increase the efficiency of our calculations by truncating the model 
space further. We used A^^ = 15 for T < 1.5 6 and Nr = 10 for T < 0.8 6. 

4.2. Signatures of pairing correlations 

We used AFMC to calculate the heat capacity and spin susceptibility as a function 
of temperature for both even and odd number of electrons. Since the imaginary time is 
discretized, it is necessary to correct for systematic errors in the size of the time slice A(3. 
We calculated the quantities of interest for two time slices Af3 = 1/32 and Aj3 = 1/64, 
and then extrapolated the results to A/3 = 0. 

The heat capacity is calculated using the same method as in the nuclear case [6]. The 
left panel of Fig. [3] shows the heat capacity versus T for a grain with A/ 6 = 1. The 
AFMC calculations for even (odd) number of electrons are shown by solid (open) circles. 
The dashed line is the BCS result, for which the heat capacity displays a discontinuity 
at the critical temperature. In the AFMC calculations, the heat capacity is a smooth 
function of temperature, but for an even number of electrons we observe an enhancement 
of the heat capacity similar to the S'-shape heat capacity observed in even-even nuclei. 
At higher temperatures, the heat capacity approaches the Fermi gas result C = 2tt'^T/36 
valid for T ^ S (but still much below the Fermi energy). 

The spin susceptibility x = dM/dB\B=o describes the magnetic response of the grain 
(i.e., magnetization M) to a weak magnetic field B. It is given by 



gr 1 



(12) 



^ arcsinh ( ^a/V ^ ) 




(13) 



where (^s is the (^-factor of electrons in the grain, pi, is the Bohr magneton and Sz is 
the z component of the total spin of the grain. The right panel of Fig. [3] shows the 
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Figure 3. Left: the heat capacity of an ultra-small metallic grain with A/5 = 1 versus 
temperature T. Solid (open) circles are the AFMC results for a grain with an even (odd) 
number of electrons. The results are compared with the BCS approximation (dashed line) 
and the canonical independent-particle model (dotted). Right: the spin susceptibility x 
(in units of Xp) versus T for a grain with A/S = 1. Symbols and lines are as in the left 
panel. After Ref. [5]. 



AFMC spin susceptibility versus temperature for both even (solid circles) and odd (open 
circles) number of electrons. The dashed line is the BCS result and the dotted lines 
are the canonical Fermi gas results. In comparison with the canonical Fermi gas results, 
we observe that pairing correlations suppress the spin susceptibility and even-odd effects 
persist to higher temperatures. The susceptibihty for an odd number of particles diverges 
in the limit T — already in the Fermi gas picture (i.e., without a pairing interaction). 
This is a result of the finite spin of the odd particle occupying the orbital at the Fermi 
energy. However, in the presence of pairing correlations, the spin susceptibility for an 
odd number of electrons initially decreased with decreasing temperatures, displaying a 
minimum at a finite temperature before diverging at T ^ 0. For large temperatures 
the spin susceptibility approaches the value Xp = obtained for a Fermi gas at 

T>5. 

The above results for the heat capacity and spin susceptibility of nanoparticles demon- 
strate that pairing correlations in the fluctuation-dominated regime manifest through 
effects that depend on the particle- number parity [2]. 



5. CONCLUSION 

We discussed pairing correlation in nuclei and nanoparticles in the crossover between 
the bulk BCS limit and the fluctuation-dominated regime. Auxiliary-field Monte Carlo 
methods are used to include fluctuations beyond the BCS mean-field theory. The heat 
capacity and moment of inertia of nuclei exhibit pairing correlation that depend on the 
particle-number parity of protons and neutrons. Similar even-odd effects (in the number 
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of electrons) are observed in the heat capacity and spin susceptibihty of nanoparticles. 
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